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ABSTRACT. 


A phylogeny of Phalaenopsis was reconstructed using three chloroplast markers, matK, atpH- 


atpF, and trnD-trnE, totaling 2177 base pairs. Results of Bayesian and Maximum Parsimony analyses 
support the placement of the species of Doritis and Kingidium into a more broadly defined Phalaenopsis, 
as proposed in a revision of the genus by Christenson. While some of Christenson’s subgeneric groups 
appear to be monophyletic, several species may need to be reclassified, if a natural classification is desired. 
Taxon sampling, however, must be completed, and the identities of several plants confirmed, before defin- 
itive conclusions can be made. Another marker of approximately 1300 base pairs will be added to the 
analysis. It is expected to provide more support and to improve the resolution of the tree, especially at the 


basal nodes and in the subgenus Polychilos. 
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INTRODUCTION 


Phalaenopsis is a genus of orchids containing 
63 species that occur throughout Southeast Asia 
and the Pacific Islands. Christenson (2000) re- 
cently revised the genus and included within it 
the species that were considered Doritis and 
Kingidium in the earlier work of Sweet (1980). 
Christenson divided the genus into five subgen- 
era, Proboscidioides, Aphyllae, Parishianae, 
Phalaenopsis, and Polychilos, based on flower 
color and pigmentation patterns, lip and callus 
structure, number of pollinia, and deciduous- 
ness. Subgenus Phalaenopsis was further divid- 
ed into four sections, Phalaenopsis, Stauroglot- 
tis, Deliciosae, and Esmeralda. Subgenus Poly- 
chilos was also divided into four sections, Po- 
lychilos, Fuscatae, Amboinenses, and Zebrinae. 

Here we present a molecular phylogenetic 
analysis of the genus Phalaenopsis, including 
the species formerly in Doritis and Kingidium, 
using DNA sequences from three regions of the 
chloroplast genome: the maturase K gene 
(matK), the ATP synthase genes atpH, atpF and 
their intergenic spacer (atpHF), and the trnD— 
trnE intergenic region (t7nDE). Relationships 
between the species are inferred and the mono- 
phyly of Christenson’s subgeneric groups are 
tested. 


MATERIALS AND METHODS 
Sampling 


Sixty-five individuals are included in the cur- 
rent data set, 47 of which are Phalaenopsis spe- 
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cies (one variety and 10 replicates), and four are 
outgroup taxa from three genera (Paraphalaen- 
opsis, Sarcoglyphis, and Amesiella) shown to be 
closely related to Phalaenopsis (Whitten pers. 
comm.). Two of the three species that were pre- 
viously placed in Doritis and seven species that 
were Kingidium have been sampled. All individ- 
uals included were obtained through horticultur- 
al sources, either private or commercial growers, 
or the collections of the Marie Selby Botanical 
Gardens, the New York Botanical Garden, or the 
University of Florida. Individuals obtained 
through private or commercial growers were 
identified as they bloomed, and vouchers have 
been placed in The University of Texas Herbar- 
ium (TEX). 


DNA Extraction, PCR Amplification, and 
Sequencing 


Total DNA was extracted either from fresh or 
silica-dried tissue using the CTAB protocol of 
Doyle and Doyle (1987). Extracts were purified 
with the QIAGEN QIAEX II Suspension kit. 
Standard PCR protocols were used (Mullis & 
Faloona 1987), and dimethyl sulfoxide (DMSO) 
was added to reactions when necessary for suc- 
cessful amplification. Primers for matK and 
trnDE were taken from Whitten et al. (2000) and 
Demesure (1995), respectively. Primers for 
atpHF were designed from the complete Maize 
chloroplast genome available in GenBank 
(Maier 1995). Amplifications were visualized on 
agarose gels and purified with QIAGEN QIA- 
quick PCR purification kits. Cycle sequencing 
reactions were performed using BigDye Termi- 
nator 3.0 and visualized on an MJ BaseStation. 
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Analyses 


Sequences were assembled and edited in Se- 
quencher 4.2. Initial alignments were carried out 
using ClustalX (Thompson et al. 1997), and then 
modified in MacClade 4.0 (Madison & Madison 
2000). All three markers were combined for all 
analyses. Parsimony analyses were run in 
PAUP* 4.0 (Swofford 2002) using a heuristic 
search with ten random-addition-replicates and 
TBR branch swapping. A bootstrap analysis was 
performed to determine branch support with 
1000 replicates and 1000 trees saved from each 
replicate. 

Bayesian analyses were performed using 
MrBayes 3.0b4 (Huelsenbeck & Ronquist 
2001). Appropriate models of evolution were 
determined separately for each marker using the 
likelihood ratio test implemented in ModelTest 
3.06 (Posada & Crandall 1998). Model Test de- 
termined that the Felsenstein 81 (F81) + gamma 
distributed rate variation model was appropriate 
for trnDE and atpHE and the general time re- 
versible model (GTR) + invariant sites + gam- 
ma distributed rate variation was appropriate for 
matK. Data were partitioned in the Bayesian 
analysis, and the appropriate model was applied 
to each partition. Four MCMC chains were run 
for 3 million generations, with one tree saved 
every 100 generations. The first 3000 genera- 
tions prior to stationarity were discarded as the 
burn-in period. 


RESULTS 


The matK region used in these analyses con- 
sisted of 1301 base pairs, the trnDE region of 
493 base pairs, and the atpHF region of 383 base 
pairs after alignment. Combined, these regions 
produced a total of 2177 aligned base pairs, 179 
of which were parsimony informative. Parsi- 
mony analyses resulted in 100,000 most-parsi- 
monious trees with a score of 470, consistency 
index (CI) of 0.679, and retention index (RI) of 
0.823. 

Results of the Bayesian analysis are summa- 
rized as a majority rule consensus tree of 27,000 
trees sampled after stationarity (FIGURE 1). Par- 
simony bootstrap values and Bayesian posterior 
probabilities are indicated above the branches. 
Individuals that have not yet bloomed (and 
therefore not positively identified) are marked 
with an asterisk (FIGURE 1). The Bayesian ma- 
jority-rule tree was very similar to the majority 
rule consensus tree from the parsimony analysis. 

The Bayesian analysis recovers two major sis- 
ter clades within Phalaenopsis. Subgenus Poly- 
chilos is represented by Clade 1, with strong 
Bayesian support and moderate bootstrap sup- 
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port (100 and 62, respectively). The second ma- 
jor clade contains a polytomy of three groups, 
Clades 2, 3, and 4. It has, however, very low 
Bayesian support (62) and bootstrap support be- 
low 50%. This group collapses in the maximum 
parsimony majority rule consensus tree, with 
Clade 2 falling into a polytomy with Clade 1 
and a third group containing Clades 3 and 4. 
This group is only recovered in 61% of the 
100,000 best trees; and in the bootstrap consen- 
sus, it collapses, resulting in a polytomy of 
Clades 1, 2, 3, and 4. 

Clade 2 contains the species in subgenera 
Aphyllae (previously Kingidium) and Probosci- 
dioides. It has very strong support with a pos- 
terior probability of 100 and a bootstrap value 
of 99. Clade 3 contains species in subgenus Pha- 
laenopsis sections Esmeralda (previously Dori- 
tis) and Deliciosae (Kingidium) and subgenus 
Parishianae (as well as P. minus, which was 
previously Kingidium). This clade has moderate 
to low support with a posterior probability of 93 
and bootstrap support of 53. Clade 4 contains 
the species in subgenus Phalaenopsis sections 
Phalaenopsis and Stauroglottis. This clade is 
also well supported with a posterior probability 
of 100 and a bootstrap value of 85. 


DISCUSSION 
Doritis and Kingidium 


Christenson’s treatment of Phalaenopsis as a 
more broadly defined genus that includes species 
of the genus Doritis and Kingidium is supported 
by the chloroplast phylogeny presented here. 
Subgenus Aphyllae and the other species that 
were previously placed in Kingidium and sub- 
genus Phalaenopsis section Esmeralda, which 
was Doritis, are well supported by both the 
Bayesian analysis (100) and the parsimony anal- 
ysis (87) to be included within Phalaenopsis. 
Although they may have morphological differ- 
ences that distinguish them from the traditional 
Phalaenopsis, these differences are apparently 
adaptations to a terrestrial habit, in the case of 
Doritis, and to colder, drier climates associated 
with high altitudes in the case of Kingidium. 


Subgenus Polychilos 


Results strongly support the monophyly of 
subgenus Polychilos, although sectional rela- 
tionships within the subgenus are still unclear 
because of a lack of resolution and low support 
for some clades. Branch lengths within the 
group are extremely short (data not shown), and 
more variable markers will be needed to discern 
the relationship between these species. Section 
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Phylogenetic tree of Phalaenopsis reconstructed with Bayesian analysis. Numbers above branches 


indicate parsimony bootstrap value/posterior probability. Taxa for which species identity has not been confirmed 


are marked with an asterisk (*). 

Amboinenses, the most species-rich section of 
subgenus Polychilos, is well sampled. There is 
support for a close relationship between Pha- 
laenopsis hieroglyphica, P. pulchra, and P. 
lueddemanniana (posterior probability of 99), 
which is not surprising, since P. hieroglyphica 
and P. pulchra were originally described as va- 
rieties of P. lueddemanniana. Phalaenopsis 
mariae and P. pallens are sister taxa with P. 


| 
| 


bastianii occurring basally within the group. 
Sweet (1980) pointed out that P. mariae and P. 
pallens have similar lip morphology, and Chris- 
tenson (2000) noted the probable, close relation- 
ship between P. mariae and the relatively re- 
cently described P. bastianii. 

Unfortunately, taxon sampling in the sections 
other than Amboinenses is not thorough. Section 
Polychilos is represented by two of the four spe- 
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cies in that group, Phalaenopsis mannii and P. 
cornu-cervi, which occur basal to most of the 
species in Amboinenses. Sections Fuscatae 1s 
represented by only one of the four species, P. 
fuscata, which occurs in a clade with P. dow- 
eryensis, P. gigantea, and P. maculata that is 
basal to the rest of subgenus Polychilos. Section 
Zebrinae is represented by two of its five spe- 
cies, P. tetraspis and P. inscriptiosenensis, 
which are well nested with the species in section 
Amboinenses. 


Subgenera Aphyllae and Proboscidioides 


Most of the species in subgenus Aphyllae oc- 
cur as a monophyletic group (Clade 2), except 
for Phalaenopsis minus. In this analysis, P. mi- 
nus occurs, with moderate support (posterior 
probability 93, bootstrap 53), basal in a group 
composed of members of subgenus Parishianae 
and subgenus Phalaenopsis sections Esmeralda 
and Deliciosae (Clade 3). In some trees it falls 
out in a paraphyletic grade between Clades 3 
and 4 (basal to Clade 4). Placement by Chris- 
tenson (2000) of P. minus in subgenus Aphyllae 
was tenuous, and he noted that it “‘is an oddball 
of sorts and does not appear to be closely related 
to the other species in the subgenus.” The po- 
sition of Phalaenopsis lowii, of the monotypic 
subgenus Proboscidioides, sister to subgenus 
Aphyllae, is strongly supported in Bayesian and 
bootstrap analysis (100 and 99, respectively). 


Subgenus Parishianae 


Subgenus Parishianae is represented by two 
of the four species in the group, Phalaenopsis 
lobbii and P. gibbosa. It was thought that three 
species were included, but an individual tenta- 
tively identified as Phalaenopsis parishii was 
identified as P. lobbii. The fourth species of sub- 
genus Parishianae, P. appendiculata, is ex- 
tremely rare in cultivation and has not been 
available for sampling. Results strongly support 
P. lobbii as sister to subgenus Phalaenopsis sec- 
tion Deliciosae and closely related to section Es- 
meralda. Although it occurs in the same clade, 
these analyses do not support P. gibbosa as most 
closely related to P. lobbii; however, the se- 
quences representing P. gibbosa are from an in- 
dividual that has not yet flowered and may be 
misidentified. 


Subgenus Phalaenopsis 


Subgenus Phalaenopsis sections Phalaenop- 
sis and Stauroglottis form a_ well-supported 
monophyletic group and perhaps should be com- 
bined to form a subgenus with no sectional de- 
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limitations. The species in section Stauroglottis 
(Phalaenopsis equestris, P. lindenii, and P. ce- 
lebensis) do not form a monophyletic group; 
rather they occur as a paraphyletic grade (in the 
majority rule consensus tree) between the basal 
species of section Phalaenopsis and the more 
derived species. Support at this branch, however, 
is low (posterior probability 82 and bootstrap 
65), and more data may recover different rela- 
tionships. 

Sections Esmeralda and Deliciosae do not ap- 
pear to belong in subgenus Phalaenopsis. They 
occur in a clade with subgenus Parishianae and 
Phalaenopsis minus. The relationship of this 
clade to the other clades in the genus is unclear, 
but it appears to be more closely related to sub- 
genera Phalaenopsis and Aphyllae than to Po- 
lychilos. Support at these nodes is very weak, 
and more data are necessary to make conclu- 
sions as to the relative relationships between any 
of the clades. 


CONCLUSIONS 


Phalaenopsis is a phylogenetic challenge. The 
preliminary results presented here provide an 
initial understanding of the evolutionary rela- 
tionships in Phalaenopsis. Phalaenopsis seems 
to have diversified relatively recently, especially 
within the subgenus Polychilos. Genetic varia- 
tion is low, especially in the chloroplast. To im- 
prove resolution and branch support, an addi- 
tional chloroplast marker, petD, consisting of ap- 
proximately 1400 base pairs, is being added to 
the data set. Inclusion of a nuclear DNA marker 
is also desirable. The internal transcribed spacer 
(ITS) of nuclear ribosomal DNA is the most 
commonly used nuclear marker for reconstruct- 
ing infrageneric phylogeny. It is problematic, 
however, in Phalaenopsis, as individuals exhibit 
multiple, paralogous types that may be phylo- 
genetically misleading or may in some cases in- 
dicate hybridization between species. In our 
search for nuclear markers, we have found sev- 
eral other regions of the nuclear genome that 
also occur in multiple copies. Additional plants 
are being sought to complete taxon sampling of 
Phalaenopsis as thoroughly as possible. 


LITERATURE CITED 


Christenson, E.A. 2000. Phalaenopsis: A Monograph. 
Timber Press, Portland, Oregon. 

Demesure, B., N. Sodzi, and R.J. Petit. 1995. A set of 
universal primers for amplification of polymor- 
phic non-coding regions of mitochondrial and 
chloroplast DNA in plants. Molec. Ecol. 4: 129— 
131. 

Doyle, J.J. and J.L. Doyle. 1987. A rapid DNA isola- 


SECOND IOCC PROCEEDINGS 27 


tion procedure for small quantities of fresh leaf 
tissue. Phytochem. Bull. 19: 11—15. 

Huelsenbeck, J.P. and E Ronquist. 2001. MRBAYES: 
Bayesian inference of phylogenetic trees. Bioin- 
formatics 17: 754-755. 

Maddison, W.P. and D.R. Maddison. 2000. MacClade 
4.0. Sinauer, Sunderland, Massachusetts. 

Maier, R.M., K. Neckermann, G.L. Igloi, and H. Kos- 
sel. 1995. Complete sequencing of the maize chlo- 
roplast genome: gene content, hotspots of diver- 
gence and fine tuning of genetic information by 
transcript editing. J. Molec. Biol. 251(5): 614— 
628. 

Mullis, K.B. and EA. Faloona. 1987. Specific synthesis 
of DNA in vitro via polymerase-catalyzed chain 
reaction. Meth. Enzymol. 155: 335-350. 

Posada, D. and K.A. Crandall. 1998. MODELTEST: 


testing the model of DNA substitution. Bioinfor- 
matics 14: 817-818. 

Sweet, H.R. 1980. The Genus Phalaenopsis. Day 
Printing Corp., Pomona, California. 

Swofford, D.L. 2002. PAUP* 4.0b10. Phylogenetic 
analyses using parsimony (*and other methods). 
Sinnauer Associates, Sunderland, Massachusetts. 

Thompson, J.D., T.J. Gibson, E Plewniak, E Jean- 
mougin, and D.G. Higgins. 1997. The ClustalX 
windows interface: flexible strategies for multiple 
alignments aided by quality analysis tools. Nucl. 
Acids Res. 25: 4876—4882. 

Whitten, W.M., N.H. Williams, and M.W. Chase. 2000. 
Subtribal and generic relationships of Maxillarieae 
(Orchidaceae) with emphasis on Stanhopeinae: 
combined molecular evidence. Amer. J. Bot. 87: 
1842—1856. 


